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I propose quantum versions of the Ziv-Zakai bounds as alternatives to the widely used quantum 
Cramer-Rao bounds for quantum parameter estimation. From a simple form of the proposed bounds, 
I derive both a "Heisenberg" error limit that scales with the average energy and a limit similar to 
the quantum Cramer-Rao bound that scales with the energy variance. These results are further 
illustrated by applying the bound to a few examples of optical phase estimation, which show that a 
quantum Ziv-Zakai bound can be much higher and thus tighter than a quantum Cramer-Rao bound 
for states with highly non-Gaussian photon-number statistics in certain regimes and also stay close 
to the latter where the latter is expected to be tight. 
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In statistics, one often has to resort to analytic bounds 
on the error to assess the performance of a parameter 
estimation technique. For the mean-square error crite- 
rion, the Cramer-Rao bounds (CRBs) are the most well 
known Although the CRBs are asymptotically tight 
in the limit of infinitely many trials, it is well known that 
the bounds can grossly underestimate the achievable er- 
ror when the likelihood function is highly non-Gaussian 
and the number of trials is limited [l|, [2[ . For such sit- 
uations, the Ziv-Zakai bounds (ZZBs), which relate the 
mean-square error to the error probability in a binary hy- 
pothesis testing problem, have been found to be superior 
alternatives in many cases 0, 0] . These bounds are of- 
ten much tighter in the highly non-Gaussian regime and 
can also follow the CRBs closely for large numbers of tri- 
als 0]- In physics, the ZZBs have also been applied to 
gravitational- wave astronomy ^] . 

The CRBs can be generalized for quantum parame- 
ter estimation, where one estimates an unknown param- 
eter such as phase shift, mirror position, time, or mag- 
netic field by measuring a quantum system such as an 

ftical beam, an atomic clock, or a spin ensemble 
Given a quantum state to be measured, the quan- 
tum CRBs (QCRBs) give error bounds that hold for 
any measurement, but since they are always less tight 
to the error than the corresponding classical CRBs 
the QCRBs share all the shortcomings of their classical 
counterparts. This is an outstanding problem in quan- 
tum metrology, as there have been many claims based 
on the QCRBs or other similarly rudimentary arguments 
about the parameter-estimation capabilities of certain 
exotic quantum states (lol - [l^ , but such claims cannot be 
justified if the bounds are not tight. Similar to the classi- 
cal case, one expects the QCRBs to be tight when many 
copies of the quantum object are available the ques- 
tion is how many. For example, Braunstein et al. found 



numerically that the CRB for phase estimation using the 
quantum state proposed by is ti ght only when the 
number of copies exceeds a threshold [13] , while Genoni 
et al. found experimentally that the QCRB for a phase- 
diffused coherent state is tight only after ^ 100 copies 
have been measured 
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In this Letter, I propose quantum Ziv-Zakai bounds 
(QZZBs) as alternatives to the QCRBs for quantum pa- 
rameter estimation. The QZZBs relate the mean-square 
error in a quantum parameter estimation problem to the 
error probability in a quantum hypothesis testing prob- 
lem, and should be contrasted with previous studies that 
consider quantum interferometry as a binary decision 
problem only 15[. To demonstrate the versatility of the 
proposed bounds, I show that a simple form of the bounds 
can produce both a "Heisenberg" error limit (H limit |1Q] ) 
that scales with the average energy and another 

limit similar to the QCRB that scales with the energy 
variance. I then illustrate these results by applying the 
bound to a few examples of optical phase estimation. An 
especially illuminating example is the state proposed by 
Rivas and Luis, the QCRB of which can be arbitrarily 
low [l2]- I show that a QZZB can be used to rule out 
any actual error scaling that is better than the H-limit 
scaling for multiple copies of this state. Beyond a certain 
number of copies, the QZZB starts to follow the QCRB 
closely, thus revealing the regime where the QCRB must 
be overly optimistic and indicating more precisely the 
asymptotic regime where the QCRB is tight. Although 
the QZZBs are also lower error bounds and not guaran- 
teed to be tight either, the study here and the usefulness 
of their classical counterparts suggest that they should 
be similarly useful for quantum parameter estimation in 
general, whenever one is suspicious about the tightness 
of the QCRBs. 

Let X be the unknown parameter, Y be the observa- 
tion, and X{Y) be an estimate of X as a function of the 
observation Y . Generalization to multiple parameters is 
possible ^ but outside the scope of this Letter. The 
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mean-square estimation error is 



dxdyPx,Y{x, v) X(y) - x 



(1) 



where Px,Yix,y) ~ PY\x{y\x)Pxix) is the joint prob- 
abihty density of X and Y, PY\x(y\x) is the observa- 
tion probabihty density, also called the likelihood func- 
tion when viewed as a function of x, and Pxix) is the 
prior probability density. A classical ZZB is given by 2] 

E > - / drrV / dx [Px{x) +Px{x + t)] 

^ Jo J-oo 

XPTe{x,X+T), (2) 

where Pre (a;, a; -I- r) is the minimum error probability of 
the binary hypothesis testing problem with hypotheses 
"Ho : X = X and Jii : X = x + observation densities 
Pyivlno) = PY\x{y\x), and Py(y|Hi) = Py|x(y|x + r), 
and prior probabilities Pq = Pr('Ho) = Px{x)/[Px{x) + 
Px{x + r)] and Pi = Pr('Hi) = 1 - Pq- V denotes the 
optional "valley-filling" operation Vf{T) = max^>o + 
v) 0: which makes the bound tighter but more difficult 
to calculate. Another version of the ZZB is 

pOO POO 

S > - / drrV / dx2min[Px{x), Px{x + t)] 

XVTf{x,X+T), (3) 

where Prg'(x, x + t) is the minimum error probability of 
the same hypothesis testing problem as before, except 
that the hypotheses are now equally likely with Pq = 
Pi 1/2. If the prior distribution Pxix) is a uniform 
window, the two bounds are equivalent 2]. For reference, 
Ref. j20jl includesproofs of these bounds, following closely 
the ones in Ref. [3- 

To apply the bounds to the quantum parameter es- 
timation problem, let px be the quantum state that 
depends on the unknown parameter X and E{Y) be 
the positive operator-valued measure (POVM) that mod- 
els the measurement. The observation density becomes 
PY\x{y\x) = iT[E{y)px\. The hypothesis testing prob- 
lem then becomes a state discrimination problem with 
the two possible states given by px and Px+r- The error 
probability is bounded by a lower limit first derived by 
Helstrom [i,[2l|: 

Pre(x,a: + T)> i(l-||Pop,. -Pip,+,||i), (4) 



where ||A||i = tr V A'^A is the trace norm. Since all 
the quantities in the integral in Eq. ([2]) are nonnega- 
tive, a lower quantum bound on the classical bound can 
be obtained by replacing FTe{x,x + r) in Eq. ([2]) with 
the right-hand side of Eq. resulting in a QZZB. For 

Plfix,X + T), 

PTf{x,x + r)>^ fl-l\\p,-p,+,\\,\ (5) 



1 

> - 
- 2 



1- v/i-P(p.,/ 



(6) 



where F is the quantum fidelity defined as F{px, Px+t) = 
(tr \/ \/lh:Px+r4Pxf ■ The inequality in Eq. (|6]) is 
proved in Ref. 21] and becomes an equality when px 
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(2) 

Px 



F 



is pure. For a product state p)^ 
nj^i F{px \ Pxl^), and Eq. © is especially convenient. 
Equations ([3]), ([5]), and ([6]) form another QZZB, which is 
much more tractable and shall be used in the remainder 
of the Letter. 

Similar to the Bayesian version of the QCRB ;7l il, 
the QZZBs allow one to compute lower quantum limits 
that hold for any measurement and estimation method 
by considering only the quantum state px and the prior 
distribution Px{x). There are, however, at least three 
significant differences between the two families of bounds: 
(1) The QZZBs are not expected to be saturable exactly 
in general, unlike the QCRBs in special cases [22], as 
the QZZBs are derived from the classical ZZBs, which 
are also not saturable usually, and the Helstrom bounds, 
which cannot be saturated for all x and t using one 
POVM. (2) While the QCRBs depend only on the in- 
finitesimal distance between px and its neighborhood 
the QZZBs depend on the distance between px and 
Px+T for all relevant values of x and r. (3) The QCRBs 
are ill-defined if px and Px {x) are not differentiable with 
respect to x, whereas the QZZBs have no such problem. 

Assume now that px is generated by the unitary evo- 
lution 



Px = exp(— i_ffX)/9exp(iiJX), 



(7) 



where H is a Hamiltonian operator and p is the ini- 
tial state. It can be shown that F{px,Px+t) ^ 
KV'I exp(— zifr)|-0) I , where M is a purification of p with 
the same energy statistics [23|. Write 1-0) in the energy 
basis as IV') = Efe Ck\Ek) with H\Ek} - Ek\Ek). Then 



F{Px,Px+t) > 



^|Cfe|2exp(-iPfeT) 

k 

ICfePlGp cos[{Ek ~ E,)t] = P(r), 



(8) 

which is independent of x. Assume further that the prior 
distribution is a uniform window with mean p, and width 
W given by 



Px (x) = — rect 



(9) 



With the optional V omitted, Eqs. dS]), dH), ([8]), and 
give 



E > = - 
- 2 



w 



(10) 



This inequality can be used to derive both an H limit and 
a QCRB-like variance-dependent limit. 
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Applying the inequality cos 61 > 1 — A|6'| to Eq. ([8]), 
where A « 0.7246 is the implicit solution of A = simp — 
(1 — cos (j))/(f> for < (/) < TT, one obtains F{t) > 



|Cfc|2|C/|2 (1 - X\Ek - Ei\t). Let Eo be the mini- 
mum Ek- Then AEk = Ek — Eg is nonnegative and 
\Ek -Ei\ = \ AEk - AEi\ < AEk + AEi, which leads to 



F(t) > 1 - 2XH+T, 



{ij\H\iP) - Eo. (11) 



A tighter bound in terms of iJ+ may be found using the 
formalism in Ref. [23] but Eq. (fTT|) suffices here. Since 
the bound in Eq. (fTTI) goes negative for r > l/(2Ai?+), 
one can use the tighter bound F{t) > there. Assuming 
a large enough H+ so that W > 1/{2\H+), Eq. (fTO]) 
becomes 



1 
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2 7o 



1/{2\H+) 
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^ for > ^ 



2Ai?4 



(12) 



Equation (|12l) is an H limit that scales with the average 
energy relative to the ground state and does not depend 
on the prior W for large H+. This result is subtly dif- 
ferent from the one in Ref. [iSl, which does not average 
the mean-square error over a prior distribution and uses 
a different method to prove the limit. The limit derived 
in Ref. p^ . on the other hand, does include prior in- 
formation and is tighter than Eq. p^ . but makes the 
additional assumptions that H has integer eigenvalues 
and W = 2tt. The H limit derived here also does not 
contradict with Ref. jlJ], which assumes H oi , fc an 
integer, and defines a different H limit in terms of n. 

To derive another limit in terms of the energy variance, 
note that the fidelity can also be bounded by [1^ [2^ 



F{t) > cos^AHt) for < r < 

AH' = {^\H^\^}-{i;\HW- 
With W > Tr/2AH, Eq. ^ becomes 



(13) 



1 



TT/2AH 



drr ( 1 - — j [1 - sin(Ai7r) 



2 7o 

7rVl6-l/2 l + 7rV48-7r/2 
AlP WAIP 

7rVl6- 1/2 „ „. TT 
AlP 



for W > 



2AH' 



(14) 



which is less tight than the QCRB Ec = 1/{4:AH^) by a 
constant factor but shows that the QZZB is also capable 
of predicting the same scaling with the energy variance. 

Consider now the problem of phase estimation using a 
harmonic oscillator, assumed here to be an optical mode, 
with H = n, the photon-number operator. For compar- 
ison, the Bayesian QCRB that includes a prior Fisher 



information 11 = / dxPx{x)\d\x\P{x) / dx]"^ is 

1 



E > E 



c 



4A7V2 + n ' 



(15) 



where AiV^ = {'ip\n'^\ip) - iV^ and N = n is 

ill-defined for the prior given by Eq. ([9]); I shall instead 
use a Gaussian prior distribution with variance W'^/12 
for the QCRB, so that H = 12/W^. For large AN^, 
the prior information is irrelevant to the QCRB. In this 
regime, Ref. [20] shows that the QZZB is less tight than 
the QCRB by just a factor of 2 when the photon-number 
distribution |CmP can be approximated as continuous 
and Gaussian, a case in which the QCRB is known to be 
saturable 22]. Thus the two bounds can differ substan- 
tially only when |Cmp is highly non-Gaussian. 

Consider first a coherent state = 

exp(-iV/2)X;™=o(^'"^V\^)l"^> witli mean pho- 
ton number N. AN'^ — N, and the fidelity is 
F{t,N) = exp[2iV(cosr — 1)], as shown in Fig. [Ija) for 
some different iV's. For a product of coherent states, 
nj=i^(T,^j) = ^(T,Ej=i^i) is identical to that for 
one coherent state with the same total photon number 
on average. 

For W = 2tt, it can be shown \2d\ that Eq. ^ gives 



1-3/2 



8VN 



exp(-47V)erfi(2V^), (16) 



where erfiz = {2/y/n) duexp{u'^). The QZZB and the 
QCRB are plotted in Fig.IUJb). In the limit of iV > 1, the 
right-hand side of Eq. ([T6|) approaches 7r/(16A^), which 
is slightly less than the QCRB given by Ec l/(4iV) 
but still obeys the expected 1/N "shot-noise" scaling. 

Next, consider the state {i/j) = {M + 1)^1/2 ^^^^^ jm), 
which has an equal superposition of number states up 
to \M) and shall be called the rectangle state here, with 
N = M/2 and AiV^ = M{M + 2)/12. The QCRB given 
by 



Sf7 = 



1 



47V(A^+l)/3 + n 



(17) 



follows the H-limit scaling 1/A^^ for large N. The fidelity 
isF(r) =sin2[(2A^-hl)T/2]/[(2A^+l)2sin2(T/2)]. Unlike 
the coherent states, the fidelities for the rectangle states 
have sidelobes, as shown in Fig. [Ijc). 
The QZZB for = 27r is [20] 



Ez > E', 



27V 

TT X ^ 1 

^ ^ 2(27V+ 1)2 2fc-H 1 
7rln(4iV-|- 1) 



(18) 



approaching a slower hiN/N"^ scaling for large N. The 
additional factor of In makes the QZZB diverge from 
the QCRB, as shown in Fig.IlJd). The InN/N"^ scahng 
was also observed previously for the phase-squeezed state 
using other methods [26j . 
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(a) Coherent States (JV = 0. 1 10) 



(b) Lower Bounds for Coherent States 




FIG. 1: (Color online). Left column: fidelity of a pair of 
(a) coherent states, (c) rectangle states, and (e) products of 
Rivas-Luis states with e = 0.1 and Nj = 1 with phase differ- 
ence T. Right column: mean-square-error lower bounds versus 
the average photon number on a log-log scale for (b) coherent 
states, (d) rectangle states, and (f) products of Rivas-Luis 
states with e = 0.1, A'j — 1, and W = 2tt; straight lines 
connecting the numerically calculated points are guides for 
eyes. 



As the final example, consider the superposition of the 
vacuum with a state |C) that has a large photon- number 
variance, viz., 1-0)^ = co|0) -I- ci\() with |cip <C 1, as 



proposed by Rivas and Luis [IJ 



can be rewritten 



as ~ Vl — e|0) + y/e\^s), where jV's) is |C) minus the 
vacuum component and renormalized. If |V's) has a mean 
photon number and photon-number variance given by 
"fN^ with 7 a constant, the mean and variance for 
are A^^- = eN^ and A7V| = [(1 + j)/e~ l]Nf. ANf can 
be made arbitrarily larger than Nj by reducing e, and 
the QCRB can be made arbitrarily small. With i' copies 
and a total photon number TV = vNj , 



1 



can decrease faster than the H- limit scaling if ve 

also decreases with N [12i] . 

The fidelity tells a very different story. Defining the 
fidelity for lips) as Fg, the fidelity for is Fj = (1 — 
e)2 + 2e(l-e)(V's|cosnT|i/'s) + e^i^s > (l-e)2-2e(l-e) = 
1 — 4e -I- 3e^. Regardless of Fj is bounded from 

below by a constant close to 1 if e ^ 1. For i/ copies, 
F = F^" > {1 - Ae + Se'^y, and a bound on the QZZB 
follows: 



12 



1 - (1 - 4e + 3e2 



(20) 



'i[{l + j)/{i^e)-l/iy]N^ 



n 



(19) 



This bound means that the actual error cannot deviate 
substantially from the prior value VF^/12 until i^e ^ 1, by 
which point even if the error catches up with the QCRB 
it can no longer beat the 1/iV^ scaling. This result is 
unsurprising in light of the now proven H limit. 

To study the behavior of the Rivas-Luis state in more 

detail, let - VT~e\0) + y^E™=i 1"^)- Fig- 

ure[IJe) plots the fidelities for some products of the Rivas- 
Luis states with e = 0.1 and iV = 1, showing sharp fea- 
tures due to lips) near r = but quickly dropping off 
to the nonzero backgrounds due to |0). Figure [Hf) plots 
the QZZB (calculated by numerically integrating Eq. (fTO|) 
with W = 2tt) and the QCRB given by Eq. versus 
the total photon number N. The QZZB is much higher 
than the QCRB for small N and comes down only when 
iV > 10 and i^e > 1. The QZZB then reaches a threshold, 
beyond which it follows closely the QCRB. This threshold 
behavior is encountered frequently in classical parameter 
estimation [l], Q and also observed in a numerical study 
of quantum phase estimation [l^ . 

In conclusion, the QZZBs are shown to be versatile 
error bounds that can predict different types of quan- 
tum limits using one unified formalism and can be much 
tighter than the popular QCRB for optical phase esti- 
mation in certain cases. To model quantum sensors more 
realistically, the QZZBs may also be generalized for wave- 
form estimation in a way similar to the QCRB @, if an 
error bound for continuous quantum hypothesis testing 
(27j can be found. 

Discussions with L. Maccone, V. Giovannetti, 
M.J.W. Hall, S. Lloyd, L. Davidovich, B.M. Escher, 
R.L. de Matos Filho, M.G.A. Paris, A. Rivas, A. Luis, 
S. Cuba, A. Tacla, CM. Caves, L.C. Kwek, A. Ling, and 
J. P. Dowling are gratefully acknowledged. This work is 
supported by the Singapore National Research Founda- 
tion under NRF Grant No. NRF-NRFF2011-07. 
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Appendix A: Supplementary Material 

This Supplementary Material contains detailed derivations of some of the results presented in the main text. 
Section lA II derives the classical Ziv-Zakai bounds, Sec. lA 21 calculates a quantum Ziv-Zakai bound (QZZB) when 
the photon-number distribution can be approximated as continuous and Gaussian, Sec. lA 31 calculates the bound for 
coherent states, and Sec. lA 41 calculates the bound for rectangle states. 



Let 



1. Derivation of the classical Ziv-Zakai bounds 

a^\X{Y)-X\ 



(Al) 
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be a nonnegative random variable. The mean-square error becomes 

POO 

S = / dsP„{s)s^, 
Jo 

where Pa{s) is the probabiUty density of a. With 

and integration by parts, S becomes 

E = 2 / dssPT{a>s) = - drrFT (\X - X\ > -) . 
Jo 2 Jo ^ 2 / 

Since both r and Pr(|X — X| > t/2) are nonnegative, one can find a lower bound on E by bounding Pr(|X- 
Rewrite Pr(|X -X\> t/2) as follows: 

Pr(|X-X|>I)=Pr(A--X>I)+P,(i-X<-I) 



J — c 



rfa;oPx(a;o)Pr(X-X> - 



J — c 



X = Xq 



and let 



This yields 



= X, 



+ / dxiPx(xi)Pr(X-X<-- 



a;i = a; + T. 



Pr {\X - X| > = y"°° Px{x) Pr 



Px(a;) Pr(X>a;+- 



X = x 









X = a; + 





/oo 
da;[Px(a;) + Px(a; + T)] 
-OO 



X {PqV-s:[X> X 



X = x\ + PxI>r[X <x + 



X =X + T 



where 



Po 



Px{x) 



Pi 



Px{x + t) 



Pxix)+Px{x + T)' Px{x)+Px{x + t) 

Now consider a binary hypothesis testing problem with hypotheses 



l-Po- 



Ho X = X, 

prior probabilities 

Pr(Ho) = Po, 
observation probability densities 

PY{y\no) = PYix{y\x), 

and the following suboptimal decision rule: 



Hi:X = x + T, 
Pr(^i) = Pi, 

PY{y\'Hi) = PYix{y\x + T), 



(A2) 

(A3) 

(A4) 
-X\ > t/2). 

(A5) 

(A6) 
(A7) 



(A8) 

(A9) 

(AlO) 

(All) 

(A12) 
(A13) 



choose Ho 'ifX < x + —, 



choose Hiif X > x + —. 



(A14) 
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The error probability of this hypothesis testing problem is then precisely given by the expression in the curly brackets 
in Eq. (|A9|) . This expression is lower-bounded by the minimum error probability of the hypothesis testing problem, 
denoted by Pre (a;, x + t), which does not depend on X, and one obtains 

Pr(|X-X| > > y" dx[Px{x)+Px{x + T)]PT,{x,x + T). (A15) 

The left-hand side is a monotonically decreasing function of r, so a tighter bound can be obtained if we fill the valleys 
of the right-hand side as a function of r. Denoting this valley-filling operation as V: 



one gets 



V/(t) = max/(r-f7y), (A16) 

Pr(|X-X| > ^) > Vy dx[Px{x)+Px{x + T)]Pr,{x,x + T), (A17) 
S>:r/ drrV dx[Px{x) + Px{x + t)] Pi e{x,x + t). (A18) 

This is a Ziv-Zakai bound. Another version that relates the mean-square error to an equally-likely-hypothesis-testing 
problem can be obtained from Eq. (jA8| : 



Pr (\X - X\>^) > I dx2m:in[Px{x),Px{x + T)] 

^ ~ X ^ x + t] \ . (A19) 



X < - Pi \ X > X + - 



X = x \ + \pi \ X < x+ - 



The expression in the curly brackets is now the error probability of the same hypothesis testing problem as before, 
except that the prior probabilities are 

Pr(Ho) = Pr(Hi) = i. (A20) 

Another Ziv-Zakai bound follows: 

^>7: drrV dx2mm[Px{x), Px{x + T)]PTf {x, x + t), (A21) 

2 Jo J-oo 

where Prg'(a;, a; -t- r) now denotes the minimum error probability with equally likely hypotheses. Note that these 
bounds make no assumption about the estimate. 



2. Quantum Ziv-Zakai bound for approximately Gaussian photon-number distributions 



Consider the QZZB for optical phase estimation with a uniform prior window: 



1 



w 



Ez = - / drr 1 

F{t) = |^|C„|2exp(znT) 



(A22) 
(A23) 



where |Cnp is the photon- number distribution. If |C„p can be approximated as a Gaussian distribution and the sum 
in Eq. (|A23p as a continuous Fourier transform. 



I a 



Fir) 



1 



2ttANj 



exp 



2ANf 



dn\Cn\'^ exp(mT) 



expl-AN-T 



(A24) 
(A25) 
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Since 1 - > F/2, 



dTT 1 



for W ^ oo, 



(A26) 



which is lower than the quantum Cramer-Rao bound (QCRB) by a constant factor of 2. This shows that the two 
bounds can differ significantly only when |C„p cannot be well approximated by a Gaussian. 
For multiple copies, the fidelity can be written as 



F(r) 



(A27) 



which is the squared magnitude of the Fourier transform with respect to the total photon number nj . By virtue 
of the central limit theorem, the total-photon-number statistics will become approximately Gaussian with variance 
i/ANj in the limit of large i^. This means that, regardless of the form of |C„p, the QZZB and the QCRB will become 
comparable in the limit of large v. 



3. Quantum Ziv-Zakai bound for coherent states 



With the inequalities 



, , T \ W . TTT 

r 1 > — sm — , 

1 



1 - \/]~F > -F, 
- 2 ' 



0<T<W, 

< F < 1, 



the QZZB becomes 



(A28) 
(A29) 



(A30) 
(A31) 



For a coherent state and W = 2tt, 

/•27r 

> - exp(-2iV) / dr sin — exp {2N cos r) . 
8 Jo 2 

Changing the integration variable to u = cos(r/2) and using the identity cost — 2cos^(t/2) — 1, one obtains 

> Trexp(-4iV) / dMexp(47Vw2), (A33) 
2 J-i 

which leads to Eq. (16) in the text. 



(A32) 



4. Quantum Ziv-Zakai bound for rectangle states 

Using Eq. (|A31|) with W = 27r, one obtains 



hi = 



TT 



8(Af + 1)2' 

^ sin2(M + l)T/2 

ar : T- . 

Q smT/2 



It can be shown using trigonometric identities that 



sm 



{M + 1)t _ . - l)r 



= sm 



, . T Mt . (M - l)r 

4 sm — cos sm ^ — 

2 2 2 



4 sin - cos 



Mt 



(A34) 
(A35) 

(A36) 



The integral Im then satisfies the recursive relation 



Im = Im-2 + 



16M 



(2M- 1)(2M + 1) 



= Im-2 + 4 



1 1 
+ 



2(M-1) + 1 2M + 1_ 



with 7_i = and /q = 4. Hence 



M ^ 



fe=o 



For large M, the discrete sum can be approximated by an integral: 



Jo 



1 



/m«4/ dA:^^^ =21n(2M+l). 



(A37) 
(A38) 

(A39) 

(A40) 



